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I.  INTRODUCTION 

The  research  pursued  under  Contract  No.  F  49620-80-C-0017  was  undertaken 

as  a  joint  project  by  D.R.  Herschbach  at  Harvard  and  J.W,  Duff,  M.B.  Faist,  and 

C.E.  Kolb  at  Aerodyne  Research,  Inc.*  in  Bedford,  MA.  During  the  course  of  the 

work,  Faist  moved  to  the  Radian  Corporation,  in  Austin,  TX  but  continued  to  take 

part  in  the  project.  Although  the  Contract  was  awarded  to  and  administered  by 

Harvard,  it  included  a  subcontract  to  Aerodyne  Research. 

The  motivation  for  this  project  came  from  the  increasing  interest  in 

high  energy  or  hyperthermal  (E*^0.5  to  10  eV)  inelastic  and  reactive  collisions 

over  the  past  several  years.  This  interest  has  been  stimulated  for  at  least  two 

1  2-4 

reasons.  First,  the  advent  of  seeded  supersonic  and  arc  heated  nozzle  beams 

5  *6 

as  well  as  fast  ion  sputtered  and  laser  blowoff  beam  sources  has  opened  up 
this  previously  inaccessible  experimental  domain.  These  techniques  may  be  used 
to  produce  not  only  high  energy  closed-shell  species  (e.g.,  rare  gases-supersoni c 

nozzle^)  and  monovalent  species  (e.g.,  alkali-sputtering^  and  supersonic  nozzles 

2  3 

and  H-atom  arc  heated  sources  )  but  also  polyvalent  nonmetals  (e.g.,  N-atoms 

4 

and  0-atoms  via  discharge  or  arc  heated  nozzles)  and  polyvalent  metal  atoms 
(heated  nozzle  and  laser  blowoff  sources).  Therefore,  a  vast  range  of  funda¬ 
mental  experimental  information  involving  hyperthermal  inelastic  and  reactive 
atom/molecule  scattering  is  becoming  accessible  to  chemical  physicists. 

Secondly,  high  energy  collisions  are  of  interest  from  a  practical  as 


well  as  intellectual  point  of  view.  For  instance,  powered  missile  flight  in  the 
Earth's  atmosphere  above  90  km  is  generally  characterized  by  very  high  levels 
of  infrared  radiation  emanating  from  the  rocket's  exhaust  plume.  The  bulk  of 
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this  radiation  is  due  to  infrared  active  vibrat ional/rotat Tonal  and  pdre 
rotational  modes  of  exhaust  gases  such  as  CO^ ,  H^O,  HF,  HC1 ,  CO,  and  NH^  which 
have  been  populated  well  above  local  thermodynami c  equi l i br i um  (LTE)  levels 

7.  g 

in  hyperthermal  collisions.  The  high  energy  nature  of  exhaust  species/ 
atmospheric  species  (and  less  predominantly  exhaust  species/exhaust  species) 
collisions  in  high  altitude  rocket  plumes  is  a  result  of  the  high  velocities 
(2-6  km/sec)  achieved  by  ballistic  missiles  at  high  altitudes. 

The  development  of  sophisticated,  spaceborne  optical  tracking  and 
surveillance  systems  utilizing  sensitive  infrared  sensors  has  led  to  an 
increased  interest  in  sources  of  mi ss i le- re  1 ated  IR  signature.  Recently,  pre¬ 
dictive  computer  codes  have  been  developed  at  Aerodyne  Research  to  calculate 

o 

rocket  exhaust  plume  infrared  radiation  signatures  observable  to  both  onboard' 
and  offboard  sensors.*^  The  construction  of  these  codes  has  required  compu¬ 
tational  advances  in  nonequilibrium  fluid  mechanics,1*  *^  inelastic  energy 

transfer  in  hyperthermal  collisions,***  *^  and  nonequilibrium  molecular 
18-21 

radiation.  The  most  difficult  of  these  calcul ational  efforts  is  the 

determination  of  inelastic  excitation  cross  sections  for  high  energy  (>1  eV) 
collisions  between  common  infrared  active  exhaust  species  (e.g.,  H^O  and  CO^) 
and  atomi c  .oxygen ,  which  is  the  dominant  atmospheric  species  in  most  of  the 
upper  atmosphere. 

In  view  of  the  importance  of  these  high  energy  collision  processes  it 
Is  clear  that  there  is  a  need  for  practical  theoretical  methods  to  aid  in  the 
understanding  and  prediction  of  the  outcome  of  such  collisions. 


II.  SUMMARY  OF  RESULTS 


The  part  of  the  work  conducted  at  Harvard  has  produced  three  published 
22-23-24 

papers.  These  deal  with: 

(1)  A  "diatomics-in-molecules"  approximation  for  constructing  potential 
energy  surfaces  of  excited  electronic  states  of  the  prototypical  triatomic 
hydrogen  system. 

(2)  A  new  approximation  method  to  estimate  cross  sections  for  inelastic 
or  reactive  collision  processes  from  a  "rotating  rigid  rod"  model. 

(3)  A  comprehensive  statistical  theory  for  vector  properties  of  collisions, 
Including  all  angular  distributions  describing  the  orientations  of  the  initial 
and  final  relative  velocity  vectors  and/or  the  reactant  and  product  rotational 
angular  momentum  vectors. 

The  work  at  Harvard  also  led  to  a  discovery  of  great  promise: 

(4)  A  new  and  novel  "dimensional  interpolation"  approach  to  accurate 
electronic  structure  calculations.  In  this,  calculations  are  carried  out  for 
the  rather  simple  limiting  cases  correspond! ng  to  one-dimension  (D  «  l)  and 
infinite  dimension  (D  a  «).  Interpolation  in  terms  of  1/D  then  provides  the 
electronic  energy  for  the  real  system  (D  =*  3)-  For  the  helium  atom  and  other 
two-electron  systems,  this  method  has  given  remarkably  accurate  results  (error 
less  than  5  x  10  with  little  labor.  If  this  procedure  proves  applicable  to 
larger  systems,  It  will  revolut ionize  electronic  structure  calculations. 

The  part  of  the  work  conducted  at  Aerodyne  has  produced  two  noteworthy 

results: 

(5)  A  satisfactory  theoretical  treatment  has  been  achieved  which  accounts 
for  "ballistic  energy  transfer."  This  denotes  extremely  efficient  inelastic 


transfer  of  translational  to  rovibrational  energy  observed  in  hyperthermal 

25 

collisions  of  rare  gas  atoms  with  alkali  halide  molecules.  A  very  detailed 
analysis  was  carried  out  for  the  Ar-Cs-F  system,  by  constructing  a  potential 
energy  surface  by  the  di atomi cs* i n-molecul es  method  and  computing  quasi¬ 
class  i  ca  1  col  1  i  s  i  on  t  rajec'tori  es  for  a  wide  range  of  conditions.  The  results 
are  in  excellent  agreement  with  experiment.  Thus  they  serve  to  account  for 
ballistic  energy  transfer,  a  mechanism  which  had  eluded  theoretical  explanation 
for  over  ten  years. 

(6)  A  potential  energy  surface  has  been  constructed  for  the  H  +  CO 

system  by  combining  results  from  ab-initio  calculations  with  spectroscopic  data 

on  the  HCO  radical.  Trajectory  calculations  employing  this  surface  were  carried 

out  to  compare  with  the  exceptionally  good  experimental  results  recently 

obtained  for  hyperthermal  inelastic  collisions  of  H  +  CO.  Excellent  agreement 

26 

was  found  for  the  data  on  t rans 1  at i onal - to- vi brat i onal  energy  transfer. 

Agreement  with  the  data  on  vibrational- to- rotat ional  transfer  is  also  good  in  an 

average  sense.  However,  the  calculations  failed  to  predict  a  pronounced 

27 

bimodal  structure  observed  in  the  energy  transfer,  and  this  appears  to  be  due 
to  Inaccuracy  in  an  Important  region  of  the  surface. 

Three  papers  describing  the  work  outlined  in  (l),  (2),  and  (3)  are  attached 
to  this  report.  Papers  presenting  the  work  of  (*t),  (5),  and  (6)  are  still  in 
preparation.  The  funds  provided  by  the  Contract  (and  its  subcontract)  were 
fully  expended. 
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III.  A  REMARK 

Theoretical  methods  for  practical  calculations  of  collision  properties 
are  now  available,  but  cannot  be  usefully  applied  because  of  the  lack  of 
sufficiently  accurate  potential  energy  surfaces.  Methods  to  compute  such 
accurate  surfaces  thus  deserve  the  highest  priority  if  progress  is  to  be  made 


in  this  area. 
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